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X) ■ ABSTRACT 



We present results from a careful and detailed analysis of the structural and 
dynamical properties of a sample of 29 disc-like objects identified at z = 
\ in three AP3M-SPH fully consistent cosmological simulations. These simu- 

lations are realizations of a CDM hierarchical model, where an inefficient 
[ Schmidt law-like algorithm to model the stellar formation process has been 

^ ' implemented. We focus on properties that can be constrained with available 

data from observations of spiral galaxies, namely, the bulge and disc struc- 
tural parameters and the rotation curves. Comparison with data from Broeils 
(1992), de Jong (1996) and Courteau (1996, 1997) gives satisfactory agree- 
^— ^ ' ment, in contrast with previous findings using other codes. This suggests that 

the stellar formation implementation we have used has succeeded in form- 
ing compact bulges that stabilize disc-like structures in the violent phases of 
[ their assembly, while in the quiescent phases the gas has cooled and collapsed 

pi i' according with the Fall & Efstathiou standard model of disc formation. 

O ' Key v^rords: galaxies: evolution - galaxies: formation - galaxies: structure - 



galaxies: spiral - hydrodynamics - methods: numerical - cosmology: theory ■ 



t3 

^ ' dark matter 

> 

>^ '■ 1 INTRODUCTION 

H ' 

5r ! Understanding how galaxies form and evolve into the objects we observe today remains 
one of the most fundamental quest of astrophysical research. Even if the field is still at 
its beginnings, the use of numerical approaches to study how galaxies are assembled in a 
cosmological hierarchical scenario from primordial fluctuations, seems promising. The main 
advantage of these approaches (namely hydrodynamical simulations) is that physics is in- 
troduced at a most general level, and the dynamical processes relevant to galaxy assembly 
(i.e., collapse, gas infall, interactions, mergers, instabilities...) emerge naturally, rather than 
by assumption, and can be followed in detail. Only the star forming processes need to be 
modelled. These considerations emphasize the interest of hydrodynamical simulations as an 
outstanding working tool to learn about galaxy formation and evolution from primordial 
fluctuations. For other approaches to this problem, see Manrique et al. (in preparation), 
Firmani & Avila-Reese (2000) and the references therein. 

As discs are naturally produced in the collapse of a dissipative system embedded in a 
gravitational field, learning how to get well-behaved discs (i.e., similar to those observed in 
spirals) in such simulations is one important step towards understanding galaxy assembly 
in general. According to the standard model of disc formation (Fall & Efstathiou 1980), 
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disc-like structures with observable counterparts can form provided that the diffuse gas in 
the dark matter halo cools and collapses conserving its specific angular momentum, j. In 
hierarchical scenarios, halo gas cooling occurs with j conservation in the quiescent phases of 
galaxy assembly, producing extended, thin discs with an exponential mass density profile in 
the equatorial plane (Dalcanton, Spergel & Summers 1997). But violent phases play also a 
fundamental role in these scenarios. During these phases, the system undergoes interactions 
and merger events, that can destroy discs that had possibly formed in previous quiescent 
phases. 

Thus far, no hydrodynamical simulation of galaxy formation in fully consistent hierar- 
chical cosmological scenarios had provided extended discs, with structural and dynamical 
properties similar to those observed in spiral galaxies (see, however, recent work by Sommer- 
Larsen & Dolgov 1999). The objects formed were cither a) gaseous and too concentrated, 
when no stellar formation processes are considered at all (Navarro & Benz 1991; Navarro, 
Frenk & White 1995; Evrard, Summers & Davis 1994; Vedel, Hellsten & Sommer-Larsen 
1994; Navarro & Steinmetz 1997; Sommer-Larsen, Gelato & Vedel 1999; Weil, Eke & Ef- 
stathiou 1998 and references therein), or b) rather stellar and without discs, when these 
processes have been included (Navarro & Steinmetz 2000; Thacker & Couchman 2000). In 
the former cases, objects are too concentrated due to an excessive j loss in violent events; 
this is the so-called disc angular momentum catastrophe [DAMC]. In the latter cases, as 
discussed by Tissera, Saiz & Dommguez-Tenreiro (in preparation, hereafter TSDT) the par- 
ticular star formation implementations used lead to a too early gas transformation into 
stars, leaving no gas to cool and form discs at lower z. The effects of star formation have 
also been considered by Katz (1992) and Steinmetz & Miillcr (1995) who obtained, in both 
three component system that resembles a spiral galaxy. However, their simulations 
treat only the case of a semi-cosmological modellization of the collapse of a constant den- 
sity perturbation in solid-body rotation. Sommer-Larsen & Dolgov (1999) have succeeded 
in avoiding excessive j loss in warm dark matter scenarios. Nevertheless, their discs still do 
not match observations very well. 

This inability to reproduce observed discs in fully consistent cosmological simulations, 
represents clearly a major pitfall of present-day numerical simulations as a working tool 
to learn about galaxy formation and evolution. In a series of papers (Domfnguez-Tcnreiro, 
Tissera & Saiz 1998; Saiz, Tissera & Dominguez-Tenreiro 1999; TSDT) it has been shown 
that the object properties at the end of violent phases of evolution critically depend on 
the structural stability of the disc-like objects involved in interaction and merger events, 
on the one hand, and on the availability of gas to regenerate a disc structure after the last 
violent episode, on the other hand. Concerning stability, cold thin discs are known to be 
very unstable against non-axisymmetric modes (Athanassoula & Sellwood 1986; Binney & 
Tremaine 1987; Barnes & Hernquist 1991, 1992; Martinet 1995 and references therein; Mihos 
& Hernquist 1994, 1996; Mo, Mao & White 1998). Massive haloes can stabilize discs, but not 
every halo is able to stabilize the pure exponential disc that would form from j-conserving 
diffuse gas collapse in its last quiescent phase. In these central compact bulge is 

needed to ensure stabihty (Christodoulou, Shlosman & Tohline 1995; van den Bosch 1998, 
2000). Otherwise, non-axisymmetric instabilities could be easily triggered by interactions 
and mergers, causing important j losses (that is, DAMCs), followed by strong gas inflow. 

In that series of papers we also show that compact central bulges form when star-forming 
processes are considered in hierarchical hydrodynamical simulations, provided that the par- 
ticular star formation implementation used does not lead to a too early gas exhaustion (i.e., 
inefficient star formation, see Silk 1999). The physical processes involved in disc formation 
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and stability when an inefficient Schmidt law-hke stellar formation algorithm is implemented 
in these simulations are discussed in detail in TSDT. 

The aim of this paper is to show that discs identified in galaxy-like objects produced 
in these hierarchical hydrodynamical simulations are, at a structural and dynamical level, 
similar to those observed in spirals. To this end, we analyze the structural and dynamical 
properties of a sample of 29 disc-like objects formed in three different such simulations. 
We focus on those properties that can be constrained with available data for observed 
spiral galaxies. These data are taken principally from the compilations of galaxy structural 
parameters of Broeils (1992, hereafter B92), de Jong (1996) and Courteau (1996, 1997b, 
hereafter C97b). 

The paper is organized as follows: ^ describes the simulations. The general characteris- 
tics of disc-like objects at z = are outlined in and their bulge-disc decomposition and 
rotation curves are discussed in §| and respectively. Finally, a summary and conclusions 
are presented in §H. 



2 THE SIMULATIONS 

We consider three simulations, SI, Sib, and SlcQ, from different realizations of a given 
model for the primordial spectrum of the density fluctuation field. In each case, we follow 
the evolution of 64^ gas plus dark matter particles in a periodic box of comoving side 10 Mpc 
{Ho = 100 h km s~^ Mpc^^ with h = 0.5) using a SPH code coupled to a high resolution 
AP3M code (Thomas & Couchman 1992; see Tissera, Lambas & Abadi 1997 for details 
on SPH algorithm implementation). The initial distributions of positions and velocities of 
the 64^ particles are consistent with a standard CDM cosmology, with fib = 0.1, A = 
and b = 2.5; one out of ten of them is randomly chosen to be a gas particle. These initial 
distributions are set by means of the ACTION algorithm, kindly provided by E. Bertschinger. 
All, dark, gas and stellar particles have the same mass, 2.6 x 10^ Mq. AP3M-SPH integrations 
were carried out from z = 10 to the present using fixed time steps At = 1.3 x 10'' years. 
The gravitational softening length is 3 kpc and the minimum allowed smoothing length is 
1.5 kpc. These simulations include a star formation prescription described by Tissera et 
al. (1997; see Tissera 2000 for a discussion on the resulting star formation history). The 
gas cools due to radiative cooling. We use the approximation for the cooling function by 
Dalgarno & McCray (1972) for a primordial mixture of hydrogen and helium. Gas particles 
are transformed into stars if they are cold (T* < 3 x 10^ K), dense (p > pcrit ~ 7 x 10^^^ 
g cm~^), in a convergent flow and satisfy Jean's instability criterion. When a gas particle 
satisfies these conditions, it is transformed into a star particle after a time interval (r) over 
which its gas mass is being converted into stars according to 

dpgas _ dp stars _ Pgas \ 



where c is the star formation efficiency, t^, is a characteristic time-scale assumed to be pro- 
portional to the dynamical time of the particle (t* oc PgaV^) and r = In 100 t*/c is the time 
interval over which 99 per cent of the gas mass in a particle is expected to be transformed into 
stars (Navarro & White 1993). We adopt a low star formation [SF] efficiency with c = 0.01 
(for a discussion of the inefficient SF implementation, see TSDT). No supernovae explosion 
effects have been considered. However, this low c value avoids the quick depletion of gas at 

* Simulation 81 has been studied in Dommguez-Tenreiro et al. (1998) and Saiz et al. (1999), and in Tissera & Domi'nguez- 
Tenreiro (1998) (their 1.2 simulation). 
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high z, leaving enough remanent gas to form disc-hke structures at low z. To some extent, 
this could mimic the effects of energy injection from supernovae explosions. The treatment 
of supernova feedback in numerical simulations is still an open and complex problem that 
remains to be properly solved. 

3 SIMULATED OBJECTS AT 2 = 

Galaxy-like objects are identified at their virial radius, r2oo, at z = 0. Dark matter haloes 
formed in the SI simulation have been extensively studied in Tissera & Dommguez-Tenreiro 
(1998, hereafter TDT98). These haloes, as well as those formed in Sib and Sic simulations, 
are resolved with a relatively high number of particles (see the number of dark mass particles, 
A^dark, per halo, in Table p. 

The main baryonic objects embedded in these dark haloes can be identified by applying a 
friends-of-friends (fof) algorithm to the baryonic particle distribution, with a linking length 
of 10 per cent the mean total interparticle separation in the box. Baryonic objects identified 
in the simulation at 2; = are galaxy-like objects (GLOs) that span a range of morphological 
characteristics: disc-like objects (hereafter, DLOs), spheroids and irregular objects. DLOs 
have external extended, populated discs, consisting mostly of gas, central stellar bulge-like 
concentration and, in some cases, stars in a thick disc. Spheroids are mostly composed of 
stars forming a relaxed, spheroidal configuration, with a small amount of gas forming a small 
central disc. Irregular objects do not have a defined morphology, and in most cases they are 
the end product of a recent merger. In this paper only DLOs will be analyzed. Only those 
GLOs with a total baryon number A^baryon = -^star + -^gas larger than N^'^^on = have been 
considered]^. This criterion has been chosen as a compromise: taking N^^^^n much higher 
than 150 would result into too few DLOs in our sample; taking N^'^^^n niuch lower than 150 
would result into DLOs composed of too few particles to be properly analyzed. 

Out of our 3 simulations, we are able to identify different numbers of discs (9, 13 and 
7, for a total sample of 29 objects) and spheroids or irregular objects (5, 7 and 11, not 
studied in this paper) following the previous mass-limit criterion. DLOs are labelled by the 
three coordinates (in 10 x mesh units) that identify the cell where the centre of mass of their 
host halo is found in the total box. When several haloes (hosting either DLOs, spheroids 
or irregulars) are found in a given cell, a Greek letter is added to distinguish among them. 
Some haloes host more than one baryonic object; asterisks in Table |I| identify haloes that 
include more than one baryonic object (either DLO, spheroid or irregular). In this case, the 
different components are distinguished by means of a capital letter in order of decreasing 
mass. For instance, DLO ^^1650 (in Sib) belongs to a massive halo, where the more massive 
objects #165A (spheroid) and #165B (irregular) are also found; DLOs #233A and #233B 
(in SI), on the one side, and #643/9A and #643/?B (in Sib), on the other side, share the 
same halo. 

The virial radius, r2oo, maximum limiting radius, -Rmax (see §|D, dark, star and gas particle 
number, A'dark, -^star and iVgas, respectively, within spheres of radii r2oo and -Rmax for each of 
the 29 DLOs of the sample are listed in Table |l]. 

In Fig. |I](a) we plot, for the z-th baryonic particle inside the halo hosting DLO #242, 

t For isolated GLOs, A''baryon is defined unambiguously as Afbaryon = -'^baryon' "^li^''^ ^baryon baryon number as given 

by the friends-of-friends algorithm for a given GLO. However, some GLOs have close satellites that are included in the friends- 
of-friends algorithm, even if these satellites are dynamically and structurally different entities from the GLO itself. In this case, 
the actual GLO baryon number is Afbarvon = Nf°^ — . Nf^ , with Nf^ the number of baryonic particles in the j-th 

J uaiyuii baryon Z__/i baryon' baryon ^ 

satellite entering in N7°^ 

^ baryon 
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the cosine of the angle formed by its position and velocity vectors, rj and Vi, versus r,, 
its distance to the DLO mass centre. In Fig. |l|(b) we plot versus r^, where j^^i is the 
component of the angular momentum per unit mass of the i-th particle parallel to j^^^ (disc 
specific total angular momentum). The solid line is Vcir{R)R, where Vch{R) is the circular 
velocity at projected distance R (see In these Figures, stars stand for stellar particles, 
circles for gas particles and open symbols for counter- rotating particles (i.e., with j^^i < 0) of 
any kind. These Figures show that most gas particles placed at rj ^ 30 kpc move coherently 
along circular trajectories (so that they have cos(rj, Vi) ~ 0, see Fig. |l](a)) on the equatorial 
plane (and so they verify cos(jj, j'^ij — 1), that is, they have — — Kir(-Ri)-Ri, with 
a small dispersion around this value (see Fig. |l](b); Ri ~ if the i-th particle is on the 
equatorial plane): the particles form a cold thin disc. In contrast with disc gas particles, gas 
particles placed at r^ ^ 30 kpc (hereafter, halo gas particles) do not show any order: roughly 
half of them are in counter-rotation (open circles), their cos(rj,t?j) G [—1, 1] and their |jz j| 
takes any value under the solid line. 

Concerning stars, we see that most of them are found at rj ^ 2 kpc, where they form 
a central compact relaxed structure, i.e., a bulge, with ji^ without any preferred direction 
and very low \jz,i\', in some DLOs, as Figs. |](a) and |l|(b) illustrate, some stars are placed 
at ^ 2 kpc in some orderly fashion, forming a sort of thick disc. Baryons forming other 
DLOs in our sample follow similar behaviour patterns. 

In Fig. we plot the specific total angular momentum at z = versus mass for dark 
haloes, jdh (open symbols), for the inner 83 per cent of the disc gas mass (i.e., the mass 
fraction enclosed by i?opt = 3.2i?d in a purely exponential disc, where R^ is the disc scale 
length, see jg (filled symbols), and for the stellar component of the DLOs in our sim- 
ulations, js (starred symbols). We see that jg is of the same order as jdh, so that these gas 
particles have collapsed conserving, on average, their angular momentum. Moreover, DLOs 
formed in our simulations are inside the box defined by observed spiral discs in this plot 
(Fall 1983). In contrast, js is much smaller than either jdh or jg, meaning that the stellar 
component in the central parts has formed out of gas that had lost an important amount of 
its angular momentum. 

A detailed description of the physical processes leading to DLO configuration at 2; = 
we have just described is given in Domfnguez-Tenreiro et al. (1998) and TSDT. Here these 
processes are only briefly summarized: 

The net effect of shocks and cooling in the quiescent phases of evolution on disordered 
halo gas particles, is that they force them to a coherent rotation with global specific angular 
momentum conservation. Consequently, in the quiescent phases of the evolution, gas particles 
tend to settle at the inner halo regions, and if the gravitational potential has some centre of 
axial (or central) symmetry at these regions, gas particles will move on a plane on circular 
orbits around this centre, forming a cold thin disc. Most stars in the central bulges have been 
formed out of gas particles that had been involved in an inflow event with high j loss; these 
compact central bulges create a gravitational potential that is spherically symmetric on scales 
of a few kpc, ensuring j conservation. This stabilizes the system against gas inflows in future 
violent phases of the evolution. In fact, in the next major merger event, those gas particles 
involved in the merger and with high j will be placed on an intermediate disc; in their turn, 
the stellar bulge of the smaller DLO involved in the event is eventually destroyed by strong 
tidal forces; then, incomplete orbital angular momentum loss puts most of its stars at some 
distance of the centre, where, after some relaxation, they form a kind of thick disc. Finally, 
most of the disc external gas particles are supplied by infall of particles, either belonging to 
baryonic clumps (satellites) or diffuse component, completing the DLO assembly. Because 
of the SF implementation used, gas in discs are not transformed into stars. Then, simulated 
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discs are mostly composed of gas. Observed discs in spirals are mostly stellar, but the stellar 
discs inherit the structural and dynamical characteristics of the gaseous discs out of which 
disc stars have formed, so that, for the sake of comparison, the value of the parameters 
describing those properties can be safely determined from the simulated gaseous discs. A 
proper implementation of star formation with supernovae feedback should lead to a self- 
regulated star formation that would allow the formation of both compact stellar bulges and 
thin stellar discs. 



4 BULGE-DISC DECOMPOSITION 

Observationally, the structural parameters of spirals are obtained from a fit to the shape 
of the surface brightness profile. The galaxy luminosity profile is often decomposed into a 
bulge (b) and disc (d) component, such as 

E^^si^* {R) = eI^""' (R) + eI^""' (R) . (2) 
Each individual profile can be modelled using the parametrization from Sersic (1968), 

R \ ^/"l 



E^^s^\R) = S^'s'^*(0)exp 



Rc 



for c = b, d, (3) 



and where n is a shape parameter taken to be n = 1 for the disc component (i.e., a pure 
exponential law) and that is left as a free parameter for the bulge component. The de Vau- 
couleurs profile for spheroids corresponds to a choice of n = 4. Sj,^^*^*(0) and Rc are the 
central surface brightness and scale length for the bulge and disc components. 

The same profile form has been fit to the projected mass density of baryons, S™'*^^(i?), 
for the DLOs in our sample. S™'^^^(i?) is calculated by averaging concentric rings centred 
at the centre of mass of each DLO. Projections on the plane normal to the total angular 
momentum of each DLO (disc plane) have been used. These projected densities are binning 
dependent and somewhat noisy. To circumvent these problems, the integrated projected mass 
density in concentric cylinders of radius R and mass M'^^^{R), 

M^^\R) = 27r (^Tr'^\R!)R!dR!, (4) 





has been used as fitting function, instead of the projected mass density itself. Cumulative 
histograms M'^^{Ri) are constructed by adding equal height steps at those positions Ri 
where baryon particles are placed. Fig. ^ shows these M'^^^{Ri) histograms for two DLOs. 
Simulated discs, as well as discs observed in spirals, do not have a sharp end. We have taken 
as maximum limiting radius for the simulated discs -Rmax = 4.5-Rd|5 which corresponds to 
a limiting or truncation radius in surface brightness of B = 26.5 mag arcsec"^ assuming 
that discs have -B(O) ~ 21.65 mag arcsec"^. As minimum limiting radius we have taken 

-Rmin = 1.5 kpc^ Because the total baryon number inside -Rmax for each object is known (see 
Table |l]; this is equivalent to knowing the total DLO mass), we are left with four free fitting 
parameters among n, Sb(0), i?b, Sd(0) and Rd- An updated version of the minuit software 
from CERN library has been used to make the fits. 

t A recent analysis by Pohlen, Dettmar & Liitticke (2000) suggests a shallower truncation radius, with -Rmax ^ SR^. 
§ The fitting determination of bulge/disc parameters is stable against changes in Hmin- These changes affect mainly the -Rb.cflF 
parameter values, that in the less favoured case, -Rmin = 0.0 kpc, change on average only by less than a 10 per cent; R^ 
variations, on the other hand, are on average lower than the corresponding Monte Carlo errors due to finite sampling (see 
Fig. Q). 
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To compare the size of bulges with different shape (i.e., with different bulge shape pa- 
rameter n) the R\y parameter is meaningless because its value, as determined from the fitting 
procedure, depends on n. In fact, when bulges with the same physical scale and different 
shape are considered, varies by orders of magnitude when n varies by a factor of some few 
units (see Eqs. (^) and (^) below). To circumvent this shortcoming, an effective bulge scale 
length or half light radius, -Rb.efr, is often used by observers. In our case, it will be defined 
through the projected bulge mass profile by the condition M^-*''(i?b,eff) = where Mb 

is the bulge mass (i.e., it is a half mass radius). -Rb.efr is related to R^, through: 

, f2,J%sV'"U (5) 



where 7 (m, u) and F (m) are the incomplete and complete gamma functions. For n = 1, 
Eq. d^) gives -^^^ = 1.678; for n = 4 one finds ^^^^ = (7.669)'^. For arbitrary n, Eq. (^ can 
be approximated by 

{2.0n - 0.324)". (6) 



Rh 

The results for the effective bulge scale length, Rh^es, and disc scale length, R^^, corre- 
sponding to the optimal bulge shape parameter, nopt, and, also, to n = 1 (exponential bulge 
profile) are given in Table |^, together with their corresponding per degree of freedom 
parameter. For comparison. Table ^ also gives the disc scale length and the pei' degree of 
freedom for a de Vaucouleurs profile (n = 4). Note that the disc scale parameter, R^, is, in 
most cases, nearly independent of the bulge shape parameter (either nopt, n = 1 or n = 4). 
Knowing the length and mass scale parameters for simulated bulges and discs allows us to 
determine their mass bulge-to-disc ratios, 

B r(2n + l) /i?b\^ S— (0) . . 

D 2 Vi?d/ S5^^'^'^(0)' ^ ^ 

These ratios are given in Table ^ This Table shows another interesting parameter charac- 
terizing DLO structure: the transition radius, -Rtrans, where E™'^'^'^(i?) = S™^^^(i?). 

In Fig. 1^ we have drawn the best fits to the integrated projected baryon mass density, 
M^'^^{R), for DLOs #242 and #643/5A. A solid line arrow marks -Rtransi a dashed line arrow 
marks the lower limit of the fitting interval -Rmin = 1-5 kpc. This Figure illustrates the good 
quality achieved in the fittings, even for the regions where R < Rmin, and hence the stability 
of the bulge/disc parameter determination against changes in -Rmin, as mentioned before. 

We see in Table ^ that bulges for most DLOs are more adequately modelled with an 
exponential profile {n = 1) than with a de Vaucouleurs profile (n = 4). This agrees with 
observational results for early and late-type spiral galaxies (e.g. Andredakis, Peletier & 
Balcells 1995; Courteau, de Jong & Broeils 1996 [CdJB96]; Moriondo, Giovanelli & Haynes 
1999 [MGH99]). Furthermore, Andredakis et al. (1995) and CdJB96 have found that the 
bulge shape parameter is correlated with spiral type. In particular, CdJB96 find that most 
Sb-Sc and all high surface brightness spiral galaxies of later type are best fitted with a 
n = 1 profile, most Sa-Sab galaxies are best modelled with a n = 2 bulge, and only a small 
fraction of late type spirals have bulges that follow the de Vaucouleurs law. This distribution 
of shape parameters, riopt, is most similar to that of DLOs in Table 0. Based on the grid 
of shape parameters in Andredakis et al. (1995) and CdJB96, most of our simulated DLOs 
seem to be of type T > 1 or greater (see nopt in Table H). The dependence of the observed 
n values on the morphological type suggests a difference in the formation and evolution of 
the bulge component. 
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The values of the i?b and parameters, corresponding to a double exponential fit, are 
plotted in Fig. |^ for all the DLOs in our sample; squares, triangles and circles correspond 
to DLOs formed in SI, Sib and Sic simulations, respectively. The error bars in the top left 
corner of the Figure are the typical dispersions in the i?b (vertical) and -Rd (horizontal) pa- 
rameters values, arising from finite sampling and calculated through Monte Carlo realizations 
of the DLO profiles, as given by Eqs. (|^) and (^), with A^baryon particles and n = 1. For com- 
parison, in Fig. I we also show the i?b and values given by CdJB96 and Courteau (1997a) 
from one-dimensional luminosity decompositions of Sb-Sc spirals (points). We get average 
bulge/disc scale length ratios of R^/Ra = 0.09 ± 0.03 for DLOs in SI, Rb/Rd = 0.07 ± 0.03 
for DLOs in Sib and R\^/ R^ = 0.12 ± 0.03 for DLOs in Sic. The ratio for the whole sample 
of 29 DLOs is R^/Rd = 0.09 ± 0.04, which compares well with the values found by CdJB96 
in the r-band, 0.13 ± 0.07, and de Jong (1996) at K alone, 0.09 ± 0.04. MGH99 also find 
Rh/Rd = 0.11 ± 0.06. In all cases, consistency between available data and our DLO sample 
is good. Note that we miss DLOs with small values of i?b and/or R^ as compared with the 
range of values quoted by these authors. It simply reflects our lack of dynamical resolution 
to resolve smaller structures. 

In Fig. ^ we compare scale lengths obtained from one- dimensional surface brightness 
decompositions with scales obtained from projected mass density decompositions with bulge 
shape parameter n = 1. As compared with the disc, the bulge is typically less prominent 
in surface brightness than in projected mass density, since its mass-to-light ratio (Tb) is 
likely to be higher than the disc mass-to-light ratio (T^). It is not clear, a priori, if this 
fact makes the previous comparisons meaningless. To find out if this is indeed the case, one 
would have to compare the values of the i?b and R^ parameters as obtained from a fit to the 
surface brightness profile of a given DLO with those obtained from a fit to its projected mass 
density profile. Since this is not possible, in order to assess the feasibility of the previous 
comparisons, we adopt the following method: 

Assuming that both Tb and Td do not depend on i?, a given DLO with double exponential 
^mass^j^^ profile characterized by values of i?b, -Rd, S™''''''(0) and SJf''''(0), would have a 
surface brightness profile, E^'s^*(i?), given by Eqs. (@) and (|^) with the same i?b and R^ 
parameters, but with E{J^'^*(0) = ^'".y^''"'' and Sj^'^*(0) = ^^^y"^. The previous comparisons 
will make sense if fittings to Monte Carlo realizations of a given Y7^^^{R) profile and to Monte 
Carlo realizations of its surface brightness counterpart, E^^^*^*(i?), with the same number of 
particles, lead to the same values of the R\^ and -Rd parameters. This test has been performed 
for Y7^^^{R) corresponding to DLOs #242 and #545 (from SI) in the sample. Results are 
shown in Fig. |^ for = 1, 3 and 5, from where we see that consistent values of i?b and R^ 
are obtained (within errors resulting from finite DLO sampling). Therefore, the mass and 
length structural parameters are decoupled, making the comparisons valid. 

This compatibility between the simulated mass scale lengths and observed luminosity 
scale lengths in the r and K bands, suggests that the projected light density (i.e., stars) 
follows the projected mass density. This could be expected if the star formation rate in a 
disc were proportional to the baryon surface mass density (see Kennicutt 1998 and Silk 1999 
for a discussion). 

The mass B/D ratios given in Table ^ are plotted in Fig. |^ versus the corresponding 
nopt- Note that no correlation is observed between B/D and nopt- These mass B/D ratios 
are considerably higher than the corresponding luminosity B/D ratios of late-type systems. 
But it is difficult to decide at the present whether the central baryon concentrations found in 
DLOs are excessive or not, as both, disc and bulge, mass-to-light ratios, would be required 
to compare with observed luminosity B/D ratios. In fact, the only disc galaxy where the 
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masses of the disc and bulge/central component can be estimated directly is the Milky Way. 
Sommer-Larsen & Dolgov (1999) obtain B/D ~ 0.2 - 0.4, for a disc scale length of 2.5 - 
3.0 kpc. For other galaxies, some estimations give Tj ^ 1-2 (in solar units, for h = 0.65; 
Courteau & Rix 1999, hereafter CR99). Values in these ranges would make the mass B/D 
ratios obtained in the simulations roughly compatible with luminosity B/D ratios from 
observations, but more accurate measures of Tb and Td are required before a conclusion 
can be reached. It is worthwhile to note that decompositions of extended H I rotation curves 
(RCs) of some galaxies from Broeils' compilation (B92), using sub-maximal discs (Bottema 
1993, 1997; CR99), yield mass B/D ratios of about 1 (see Rhee 1996). Also, MGH99 find that 
5 out of 23 Sa-Scd galaxies for which they obtain a reliable fit to the rotation curve (RC) 
for both the bulge and the disc component, have mass B/D ratios of about 1 or larger. 
Note, however, that the majority of these decompositions are ill-constrained, as different 
combinations of bulge, disc, and halo profiles yield equally good agreement with the data 
(van Albada et al. 1985, CR99). 

Current observational uncertainties on Tb and Td preclude any consistency check for our 
derived baryon mass distributions. For more stringent constraints, one must turn instead to 
dynamical observables such as RCs. 



5 DLO ROTATION CURVES 

The shapes of spiral RCs are determined by their three dimensional total (i.e., bulge, disc 
and dark matter halo) mass distribution. 

The RCs of the simulated objects have been constructed by adding up in quadrature the 
contributions to circular rotation from the bulge, the disc (both of them formed by baryons, 
either gas or stars, see so that V^bar(^) = K?!'^) + ^^^1 the dark matter halo, 

Kim(r), 

Vlir) = Vl,ir) + VUr)- (8) 
We have adopted a softened Plummer potential (Evrard et al. 1994) 



^'W = — i = bar,dm (9) 



,2,^, _ GMi(< ry 

where Mbar(< f) and Mdm(< r) are the total baryonic and dark matter masses, respectively, 
inside a sphere or radius r. The adaptive AP3M scheme produces an effective dynamical 
gravitational softening at z = of eg = 3 kpc at the central region of the DLOs. The 
RCs for several DLOs in our sample are given in Fig. |^. Points are the tangential velocity 
component of the baryonic particles (projected on the disc plane), solid lines are Vdrir), 
dashed lines are Vdm('") and point lines Vharir). Note that the halo component, Vdm('"), gives 
rise to an almost flat contribution to the RCs beyond the very inner regions. These RCs 
fit very well the tangential velocity component of the baryonic particles up to i? ~ -Rmax, 
suggesting that these particles are in rotationally supported equilibrium within the potential 
well produced by both the dark haloes and the baryons themselves (TDT98). 

On scales encompassing the baryonic object (i.e., < -Rmax), DLOs have RCs that are 
declining in their outer parts, and some of them have central spikes produced by the bulge, 
similar to those found in some observed galaxies with mass B/D 1 or larger (for example, 
NGC 6674 and NGC 7331 in Rhee 1996, or UGC 89, UGC 1013 and UGC 1238 in MGH99). 
It is customary to try to explain the shapes of RCs as a combination of various 'fundamental 
plane' parameters such as global size (i.e., a spatial and a mass scale; Rubin et al. 1985), 
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central mass concentration (Rhee & van Albada 1996), and luminosity (Persic, Salucci & 
Stel 1996). Here, the RC shapes for simulated DLOs are parametrized as a function of disc 

scale length, R^, as spatial scale, rotation velocity at R2.2 = '2.2RJ^ V2.2 = Kir(-R2.2), as 
mass size, and central mass concentration, via the maximum or peak rotation velocity, V^^^^^, 
and the radius where this is reached, -Rpeak (C97b). Some combinations of these parameters 
have also been explored, such as the logarithmic slope, LS (Casertano & van Gorkom 1991 
[CvG91]; B92; de Blok, McGaugh & van der Hulst 1996). The LS is observationally defined 
as the slope of a straight line that fits the RC in log - log scale, from _R2.2 to -Rout (this last 
point is usually taken to be the last measured RC point, but see C97b for variations on this 
operational procedure). For simulated objects in our sample, the RCs are nearly linear in log 
- log scale for R > R2.2, so that LS does not appreciably change for i?out > -^2.2- The values 
of these parameters measured in simulated DLOs have been compared with those measured 
in RCs from: a) Sb-Sc field spirals drawn from the Courteau sample, with long-slit Ha 
spectra and r-band photometry (see C97b), and b) Sb-Sc spirals from Broeils' compilation 
(B92) of extended Hi RCs with surface photometry. We keep only galaxies with V2.2 > 150 
km s~^ (only one of Broeils' galaxies has V2.2 in the range [150, 180] km s~^). 

Apart from these parameters characterizing the shapes of RCs at scales of the baryonic 
objects, other interesting parameters at halo scales (i.e., ~ r2oo), are the velocity at the virial 

radius V200 = V;ir(r2oo) and the spin parameter, A = "^'""If 'g)^''" , with Jtot, -^tot and Mtot the 

total angular momentum, energy and mass for all particles inside r2oo- Note that both V200 
and A are not observable. 

In Table ^ we give V200, ^ and the RC parameters for our DLO sample. Note that the 
A values are within their expected range for a standard CDM scenario (Warren et al. 1992; 
Dalcanton et al. 1997; Lemson & Kauffmann 1998). In Fig. || we plot R^ versus Ar2oo- A 
correlation appears in the sense that haloes with higher spin host more extended discs. In the 
framework of semianalytical models of quiescent pure (i.e., bulgeless) disc formation with j 
conservation one approximately has Rd oc Ar2oo (see Dalcanton et al. 1997; Mo et al. 1998); 
we find an important dispersion due to the loss of angular momentum by bulge particles 
involved in DAMCs in the violent phases of the evolution. The DLO values of the parameters 
Rd and V2.2 are consistent with observations (see Tables |^ and ^ and Figs. ^ and By virtue 
of selecting DLOs with A^baryon > 150, we find high RC amplitudes mostly with V2.2 > 150 
km s~^. Our simulated DLOs are either intermediate spirals (100 km s^^ < V2.2 < 180 
km s~^, 9 objects), large bright spirals (V2.2 > 180 km s~^, -Rd > 5.25 kpc, 18 objects), or 
compact bright spirals (V2.2 > 180 km s~\ -Rd < 5.25 kpc, 2 objects). Our sample does not 
include any dwarf galaxy. As expected, V2.2 is correlated with V200, albeit with a dispersion. 
This dispersion is provided by the particularities of the evolutionary history of each DLO 
in the sample, namely, the number and characteristics of the interactions and merger events 
involved in its assembly (see TDT98 for a discussion). 

The mass distribution of real spirals and their decomposition into bulge, disc and halo 
components cannot be inferred uniquely from observed RCs, as different decompositions 
yield equally good agreement with the data (van Albada et al. 1985; Broeils & Courteau 
1997; CR99). One has to postulate different scalings (mass-to-light ratios) between the bulge 
and disc components in order to infer the mass distribution of the halo. Common constraints 
for (M/L)disc include a) the so-called 'maximum disc (or maximal light) hypothesis' (van 
Albada & Sancisi 1986), which postulates that the halo dark mass component needed to fit 
the RCs should be minimum, or, equivalently, that the contribution of the disc and bulge 

^ -R2.2 is the radius at which a purely exponential disc would reach its maximum circular velocity. 
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to the inner parts of RCs should be maximum, and b) sub-maximal discs with ("y^)^ ~ 
0.60 ± 0.10. This constraint is obtained either by matching the vertical velocity dispersion, 
scale height and scale length of a thin exponential disc which yields the maximum disc 
rotation (Bottema 1993), or, independently, showing that the Tully-Fisher relation of bright 
spirals is independent of surface brightness (CR99). 

A maximum disc is defined as (^^) = 0.85 ± 0.10, at the 95 per cent confidence 

V Kir J i?2 2 

level (Sackett 1997), corresponding to ('^™-)^ = 0.72 ± 0.17. To the sub-maximal discs 

described above corresponds a ( ™ ) '^^ ^■'^^ 0.12. One notes that this choice of sub- 
maximal disc corresponds to half the amount of luminous matter at i?2.2 compared to the 
maximum disc case. 

Let us now consider our simulations. Vium,2.2 = Vium(-R2.2) can be measured directly once 
i?d is determined for each DLO (assuming that V^um = Kiar) and the different hypotheses 
above can be tested. Results are given in Table ^ for both the optimal and double exponen- 
tial fits. Vbar,2.2 IS Correlated with V2.2 (and, also, with V200, albeit with more dispersion). In 
Fig. we plot the (-^^^ ratios versus V200 for DLOs in our sample (double exponential 

fits; no substantial changes occur taking instead the values corresponding to optimal fits). 
The shaded area shows the 95 per cent interval for maximum discs according to Sackett 
(1997). Only DLOs #126 and 7^545 (from SI) are inside this area, as they have (tt^) = 

\ ^ch- J B,2.2 

0.765 and 0.763, respectively. These DLOs are somewhat peculiar: DLO #126 has a partic- 
ularly populated disc; DLO #545 has two satellites close to it. We see that most DLOs lie 
inside the allowed interval for sub-maximal discs. The mean value for the simulated DLOs 
is (^)^ = 0.67 ± 0.07, in excellent agreement with Bottema (1993, 1997) and CR99. 

In Fig. pHI only a slight correlation appears between ( ) and V200 with an important 

dispersion, in the sense that less massive haloes have slightly lower values of (17^) 

V Veil ) R2.2 

This presumably reflects the fact that less massive haloes are somewhat more centrally 
concentrated than more massive ones (Navarro, Frenk & White 1995b; see also TDT98). 
The dispersion, on the other hand, reflects the particularities of the assembly of each DLO 
in the sample. 

To quantify the importance of dark matter (or baryons) in the inner regions of DLOs 
we have used the rcross parameter (TDT98). We define Tcross as the radius where Vbar('^) ~ 
Vdm('^)5 such that RCs are dynamically dominated by baryons for r < rcross, and dark matter 
dominated r > rcross (see Table Note that most DLOs have 1 < ^^^^ < 3 and only a few 
have rcross > -R2.2; note also that there is no correlation between rcross and i?d- In Fig. ^ 



we plot versus V2.2; some correlation appears, in the sense that more massive objects 
are less dark matter dominated at their centres. This correlation disappears when V2.2 is 
normalized to the halo global velocity, V200; indicating that most of the correlation seen in 
Fig. |lT] is provided by the halo total mass. In fact, as stated above, less massive haloes are 
known to be more centrally concentrated than more massive ones. 

In Figs. 0(a) and |l^(b), we plot LS versus the disc scale lengths -Rd and versus V2.2, 
respectively, for the DLOs and, as a comparison, the objects of CvG91 compilation in the 
same dynamical range. We see that DLOs with V2.2 > 180 km s^^ do have observational 
counterparts, i.e., their LS parameters take values that are found in observed spirals with 
the same range in V2.2, even if we miss objects with LS between, say, and —0.1. Those 
with V2.2 < 180 km s~^ have excessively low LSs, but recall that only one galaxy in CvG91 
has V2.2 G [150, 180] km s~^. Moreover, according to CvG91, LSs are correlated with the 
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morphological type. Our DLOs would be T < 5, in consistence with the range in T type 
previously found from the bulge shape parameters. 

It is worthwhile to note that the value of the LS is almost independent on how baryonic 
mass is distributed inside R2.2, particularly on how it is shared between the bulge and the 
disc (it changes on average by only 10 per cent when the bulge mass is halved conserving the 
baryon number inside R2.2)- The fact that LSs are negative may reflect in part the fact that 
dark matter haloes in standard CDM scenarios are too concentrated, leading to declining 
rotation curves (Navarro, Frenk & White 1998). In fact, the declining shape of DLO RCs 
results from the fiat halo contribution to RCs (see TDT98 for a detailed discussion) and the 
dynamic baryon dominance in the inner DLO regions. Declining RCs have been observed 
in spiral galaxies. Early examples can be found in Bosma (1978, 1981), van Moorsel (1982). 
Then, CvG91, B92 and Bosma (1999), among others, also found other examples. CvG91 
suggest that this is a common feature of compact bright spiral galaxies, and also frequent 
in large bright spirals (although see B92 who shows, contrary to CvG91, that declining RCs 
do not correlate with the disc scale length). Dubinski, Mihos & Hernquist (1999) have found 
that the declining character of RCs is a necessary condition to form long tidal tails in galaxy 
mergers as those observed in some interacting galaxy systems. 

In Fig. |1^ we plot V^^^^^ versus the Rpeak/Rd ratio for the 29 DLOs in our sample. For 
comparison, we also give the peak velocities in Hi of Sb-Sc spirals from table 4.1 in RC96 
(most of them from B92 compilation), with Rpeak/Rd ratios taken from C97b. We also plot 
the V^^^^^ and Rpeak/Rd ratios measured in the Courteau sample of optical RCs (C97b). We 
confirm the previous finding that DLOs in Table have observational counterparts, even 
if the smaller ones are somewhat too concentrated as compared with Courteau data. Note, 
however, that -Rpeak is difficult to measure in DLOs and the measurement hard to assess in 
spiral galaxies. 

6 SUMMARY AND CONCLUSIONS 

We present the results of a detailed and careful comparison between the parameters char- 
acterizing the structural and dynamical properties of a sample of 29 simulated DLOs and 
those measured in observed spiral galaxies. These DLOs have been identified in three fully 
consistent hierarchical hydrodynamical simulations, where an inefficient Schmidt law-like al- 
gorithm to model the stellar formation process has been implemented. In this paper we have 
been only concerned with disc structural and dynamical properties, because they keep the 
imprints of the dynamical and hydrodynamical processes that, together with their interplay 
with star formation, are key for galaxy formation and evolution. 

The comparison of DLOs in our sample with spiral discs has been a two step work. 
First, DLO bulge-disc decomposition has been performed. This yields the bulge effective 
scale length and shape parameters, -Rb.es and n^pt, the disc scale length, R^^, and the mass 
bulge-to-disc ratio, B/D, among other DLO parameters. After having tested the robustness 
of the fitting procedure, the following results have been obtained: 

(i) The scale lengths, -Rb,efr and -Rd, and their ratio i?b,efr/-Rd, are consistent with available 
data (CdJB96; de Jong 1996; MGH99). 

(ii) The distribution of bulge shape parameters is similar to that found by CdJB96 for 
their sample of Sa-Sc galaxies. 

(iii) The mass B/D ratios are somewhat high as compared to luminosity B/D ratios, 
but current observational uncertainties on the bulge and disc mass-to-light ratios make it 
difficult to draw any conclusion from this comparison. 
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In a second step, DLO RCs have been analyzed. DLO gas particles placed at distances 
between 2 kpc and 30 kpc to the centre move roughly along circular orbits on the equatorial 
plane. They are in centrifugal supported equilibrium within the potential well produced by 
the total mass distribution (TDT98). Rotation curve shapes have been parametrized through 
i?d (or i?2.2 = 2.2 i?d) as spatial scale, V2.2 = Kir(-R2.2) as mass size, and the maximum or 
peak velocity, V^l^'^^, and the radius where this is reached, -Rpeak- These parameters have 
been compared with those measured at: a) Sb-Sc field spirals drawn from the Courteau 
sample, with long-slit Ha spectrum and r-band photometry (C97b), and b) Sb-Sc spirals 
from Broeils' compilation (B92) of extended H I RCs with surface photometry. The main 
results follow: 

(i) In contrast to findings in other fully-consistent hydrodynamical simulations (e.g., 
Navarro & Steinmetz 2000; Thacker & Couchman 2000), DLO V2.2 velocities have been 
found to be consistent with observational data. This is a consequence of disc formation with 
j conservation. 

(ii) The i?d and V2.2 values obtained indicate that DLOs in our sample are either large 
bright spirals (18), intermediate spirals (9), or compact bright spirals (2). 

(iii) The average relative contribution of baryons to V2.2 in our sample is ("y^)^ — 

0.67 ±0.07, in very good agreement with Bottema (1993, 1997) and CR99, if we take Vbar = 
Vium- Most DLOs have been found to have sub-maximal discs; only two of them have been 
found to lie inside the 95 per cent confidence interval for maximum discs according to Sackett 
(1997). The previous agreement also shows that the amount of baryon mass that has ended 
up inside R < R2.2 is not excessive, again as a consequence of j conservation. 

(iv) The ^^^^^ parameter ratio is a measure of the relative amount of dark matter at 

^ ' '"200 ^ 

the inner DLO zones. Some correlation between these amounts and the circular velocity at 
i?2.2 has been found, in the sense that less massive objects tend to be more dark matter 
dominated in their central regions. 

(v) Concerning the parameters that give a measure of the central concentrations, the 
comparison of the K?""^ 

versus plots to B92 data confirms again that DLOs have 

observational coTintcrparts, even if the smaller DLOs are somewhat too concentrated as 
compared to Courteau data (C97b). 

(vi) The logarithmic slope, LS, measures the RC slopes for R > R2.2, where, as stated, 
in most cases dark matter is already dynamically dominant. DLOs with V2.2 > 180 km s~^ 
occupy the same zone in the LS versus 1^2.2 and versus it!2.2 plots as galaxies studied by CvG91 
in the same dynamical range; those with V2.2 < 180 km s~^ do not show the tendency of less 
massive galaxies in CvG91 to have higher LS . This may be a consequence of the strongly 
concentrated mass distribution of dark haloes in standard CDM scenarios. 

A concern is in order regarding numerical resolution. DLOs are resolved with a relatively 
low number of particles. In contrast, dark matter haloes are described with a much better 
resolution. An inappropriate low gas resolution would result in an unphysical gas heating 

that could halt the gas collapse (Navarro & Steinmetz 1997). However, some works suggest 
that it is an inadequate resolution in the dark matter halo component that may produce the 
larger undesired numerical artifacts (Steinmetz & White 1997). In fact, it appears that a 
well-resolved dark matter halo, even if the number of gas particles is lower (but more than, 
say, one hundred), gives rise to a well-represented (or well-resolved) gas density profile, 
being this point the most important for both the hydrodynamics and the tracking of the 
star formation history. Moreover, to make sure that the populated and extended discs in 
the simulations do not result from unphysical gas heating or smoothing, we have run a 



© 0000 RAS, MNRAS 000, 000-000 



14 A. Sdiz et al. 

higher resolution simulation (64^ particles in a periodic box of 5 Mpc, with cosmological 
and star formation parameters similar to those in SI, Sib and Sic; hereafter HRS). Only 
one disc with mass comparable to those in these lower resolution simulations forms. It has 
1713+1380 gas+star particles and its halo 36112 dark matter particles. Its analysis has 
shown that it is populated and extended, that its structural and dynamical characteristics 
are also compatible with observations (see tabic 1 and fig. 1 in Dominguez-Tenreiro et al. 
1998, where we show that the i?b and parameters and their ratio take a typical value, 
and the specific angular momentum has been conserved) and that the physical processes 
leading to its formation are essentially the same as those that are at work in SI, Sib or Sic. 

In conclusion, the comparison between DLOs produced in our simulations and observa- 
tional data allows us to affirm that they have counterparts in the real world. This agreement 
suggests that the process operating in Fall & Efstathiou (1980) standard model for disc 
formation (i.e., gas cooling and collapse with specific angular momentum, j, conservation) 
is also at work in the quiescent phases of DLO formation in these simulations, resulting in 
discs with exponential density profiles, as predicted by Dalcanton et al. (1997). However, 
in Dominguez-Tenreiro et al. (1998) and Saiz et al. (1999) it is shown that violent episodes 
(i.e., interaction and merger events) also occur and play an important role in DLO assembly. 
In particular, it is shown that dark matter haloes formed in these simulations are not always 
able to stabilize the pure exponential disc that would form at its central region according to 
Fall & Efstathiou and Dalcanton et al. scenarios (see also van den Bosch 2000). To provide 
the right conditions for disc regeneration to occur after the last violent episode of DLO 
assembly, a compact central bulge is needed; this will ensure the axisymmetric character of 
the gravitational potential well at scales of some kpcs at all times, avoiding excessive j losses 
in violent events. In TSDT these arguments are developed, and, moreover, we prove that 
a second condition is necessary for disc regeneration: the availabihty of gas at low z, that 
is, it is necessary that the SF algorithm implementation docs not result in a too early gas 
exhaustion. The good behaviour of our DLO sample as compared with observations suggests 
that the inefficient SF algorithm used in the simulations has met both requirements. 

The global agreement with observations we have found, at a structural and dynamical 
level, also represents an important step towards making numerical approaches more widely 
used to learn about galaxy formation and evolution in a cosmological framework, i.e., from 
primordial fluctuations. These approaches will be particularly useful if, as expected, future 
improvements in numerical methods and computers speed allow to numerically simulate the 
different aspects of galaxy assembly with an increasingly high degree of realism. 
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Table 1. General Parameters 

r200 Rn 



DLO# 


''200 


-^max 






Wgas 


JVstar 








JVstar 


126 


239. 


25.9 




2709 


248 


157 




516 


185 


157 


143 


194. 


37.5 




1482 


153 


47 




473 


110 


47 


233A 


310.* 


28.5 


( 


6027 


583 


164 


) 


506 


201 


95 


233B 


310.* 


25.4 


{ 


6027 


583 


164 


) 


400 


148 


69 


242 


316. 


33.4 




6325 


621 


281 




959 


339 


278 


333 


296. 


49.7 




5213 


509 


217 




1337 


338 


212 


531 


235. 


30.9 




2573 


275 


79 




589 


224 


79 


544 


312. 


43.4 




6170 


555 


230 




1304 


349 


215 


545 


328. 


26.2 




7131 


671 


252 




794 


381 


238 



SI 



sib 



054 


252. 


35.9 




3076 


299 


138 




611 


171 


138 


110 


254. 


30.3 




3145 


268 


174 




648 


163 


174 


116a 


303. 


73.0 




5357 


455 


270 




1806 


291 


269 


116/3 


178. 


38.1 




1029 


171 


30 




460 


138 


30 


164 


222. 


29.2 




2107 


233 


62 




496 


160 


62 


165C 


467.* 


39.1 


( 


19522 


1188 


1578 


) 


556 


113 


56 


211 


186. 


43.5 




1236 


157 


12 




468 


130 


12 


212 


261. 


49.4 




3419 


226 


237 




1171 


192 


234 


633A 


312.* 


40.9 


{ 


5890 


371 


401 


) 


1068 


188 


.325 


643a 


317. 


47.4 




6176 


447 


373 




1217 


180 


358 


643/3A 


307.* 


38.1 


{ 


5583 


518 


270 


) 


695 


217 


213 


643/3B 


307.* 


32.4 


( 


5583 


518 


270 


) 


351 


146 


55 


6437 


232. 


40.8 




2409 


234 


97 




683 


181 


97 



Sic 



106 


217. 


29.5 


1966 


196 


76 


419 


122 


76 


324 


191. 


36.6 


1322 


173 


22 


395 


127 


22 


342 


280. 


34.3 


4214 


470 


114 


770 


311 


107 


444 


267. 


33.8 


3727 


352 


99 


737 


243 


99 


454 


215. 


20.4 


1907 


184 


70 


383 


149 


70 


556 


265. 


31.0 


3605 


257 


229 


707 


178 


229 


631 


215. 


18.9 


1917 


193 


59 


322 


145 


59 



Distances arc given in kpe. An asterisk indicates that the corresponding DLO is embedded in a double or multiple halo, i.e., 
that one or more other baryonic objects (either disc- like or not) exist within r2oo- This is to be taken into account on evaluating 
the total number of halo particles (numbers in parentheses). 
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Table 2. Bulge/disc parameters 



Optimal Fit 



DLO# 


"opt 


-Rb.eff 




B/D 






126 


2.6 


0.90 


6.3 


1.83 


5.72 


4.46 


143 


3.2 


0.58 


10.9 


1.45 


1.61 


4.54 


233A 


1.3 


1.00 


6.5 


1.25 


5.52 


3.57 


233B 


0.6 


1.05 


5.6 


1.22 


9.09 


2.77 


242 


1.2 


1.25 


7.5 


1.23 


15.96 


4.19 


333 


1.9 


1.36 


12.0 


1.72 


5.30 


6.59 


531 


0.9 


1.70 


6.9 


2.02 


3.85 


5.24 


544 


2.3 


1.05 


10.6 


1.79 


2.44 


5.82 


545 


1.1 


1.26 


5.9 


1.23 


3.01 


3.78 


054 


1.6 


0.89 


8.4 


2.99 


2.94 


4.80 


110 


1.8 


0.60 


7.0 


2.54 


2.64 


3.46 


116q: 


2.3 


0.48 


16.4 


1.40 


13.97 


4.12 


116/3 


1.3 


1.26 


8.7 


1.70 


1.75 


4.97 


164 


1.3 


1.19 


7.2 


2.47 


1.70 


4.83 


165C 


1.6 


0.93 


9.0 


2.37 


1.08 


4.67 


211 


1.8 


1.35 


14.0 


1.42 


0.75 


6.47 


212 


2.3 


0.53 


11.2 


1.95 


4.93 


4.04 


633A 


1.1 


0.72 


9.0 


3.01 


4.69 


3.33 


643a 


1.4 


1.19 


10.9 


3.40 


3.14 


5.91 


643/3A 


2.8 


0.86 


9.4 


2.36 


3.18 


5.75 


643/3B 


1.4 


1.28 


7.6 


1.98 


1.46 


5.02 


6437 


1.7 


1.18 


13.2 


2.43 


2.90 


6.56 


106 


2.0 


1.06 


7.5 


2.69 


1.18 


5.51 


324 


1.2 


1.42 


8.3 


1.96 


0.80 


5.22 


342 


1.3 


1.43 


8.2 


1.45 


2.52 


5.09 


444 


1.3 


1.41 


7.9 


2.02 


3.49 


5.34 


454 


1.0 


1.23 


4.6 


3.14 


2.33 


4.26 


556 


2.1 


0.86 


7.6 


3.70 


1.64 


5.37 


631 


0.9 


1.13 


4.1 


2.36 


0.97 


3.50 



Distances are given in kpc. 







n = 1 






n 


= 4 






B/D 




-^trans 




X'/N 


1.09 


5.8 


1.48 


10.40 


3.50 


6.5 


6.35 


0.91 


10.1 


1.36 


2.30 


3.52 


11.1 


1.53 


1.08 


6.4 


1.22 


5.55 


3.43 


7.1 


5.81 


0.89 


5.7 


1.23 


9.55 


2.89 


5.8 


12.54 


1.29 


7.5 


1.21 


16.19 


4.06 


7.9 


42.04 


1.50 


11.1 


1.58 


11.40 


5.36 


12.8 


11.98 


1.69 


6.9 


2.07 


4.00 


5.39 


9.3 


19.57 


1.26 


9.7 


1.61 


9.05 


4.54 


11.3 


4.06 


1.27 


5.9 


1.22 


3.11 


3.74 


6.4 


20.24 


1.04 


8.0 


2.82 


3.65 


4.13 


9.7 


4.32 


0.78 


6.8 


2.47 


3.28 


3.13 


7.2 


3.46 


0.76 


16.3 


1.39 


14.55 


3.50 


16.5 


14.79 


1.34 


8.5 


1.64 


1.93 


4.60 


10.1 


3.10 


1.25 


6.9 


2.36 


1.99 


4.44 


9.9 


4.64 


1.06 


8.7 


2.28 


1.28 


4.15 


9.8 


1.71 


1.49 


13.0 


1.36 


1.19 


5.40 


17.0 


1.30 


0.80 


11.1 


1.91 


6.30 


3.46 


11.3 


5.49 


0.74 


9.2 


2.98 


4.33 


3.32 


9.6 


6.62 


1.28 


10.6 


3.27 


5.07 


5.30 


13.2 


19.61 


1.13 


8.5 


2.02 


10.91 


4.22 


10.0 


3.70 


1.35 


7.3 


1.87 


1.96 


4.54 


9.5 


3.45 


1.32 


11.9 


2.32 


4.33 


5.30 


17.8 


5.37 


1.20 


6.6 


2.25 


2.47 


4.22 


8.2 


1.70 


1.46 


8.2 


1.91 


0.84 


4.99 


10.1 


2.56 


1.46 


7.7 


1.38 


2.87 


4.60 


9.2 


9.95 


1.46 


7.6 


1.92 


4.22 


4.89 


10.9 


11.14 


1.24 


4.6 


3.10 


2.34 


4.20 


10.2 


6.45 


1.06 


6.9 


3.17 


6.17 


4.13 


8.2 


3.70 


1.11 


4.2 


2.43 


0.97 


3.62 


18.0 


2.78 
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Table 3. Dynamical Parameters 



n = 1 Optimal Fit 



DLO# 


V200 


A 


Across 


LS 


-^peak 


X 7-p6ak 
cir 


V2.2 


Vbar,2.2 


V2.2 


Vbar,2.2 


126 


114. 


0.039 


16.6 


-0.176 


10.6 


216. 


210. 


157. 


206. 


153. 


143 


98.5 


0.045 


8.94 


-0.172 


7.98 


156. 


148. 


83.0 


147. 


80.6 


233A 


157.* 


0.063* 


14.0 


-0.139 


9.45 


203. 


197. 


139. 


197. 


139. 


233B 


157.* 


0.063* 


12.9 


-0.085 


6.93 


179. 


176. 


125. 


176. 


125. 


242 


160. 


0.069 


19.3 


-0.168 


6.33 


279. 


257. 


186. 


257. 


186. 


333 


150. 


0.079 


16.2 


-0.108 


7.25 


250. 


221. 


143. 


218. 


139. 


531 


118. 


0.011 


12.7 


-0.192 


7.70 


217. 


206. 


140. 


206. 


140. 


544 


158. 


0.023 


17.8 


-0.104 


5.69 


254. 


228. 


153. 


225. 


148. 


545 


166. 


0.049 


18.6 


-0.159 


8.66 


278. 


270. 


204. 


269. 


204. 


054 


126. 


0.122 


16.4 


-0.152 


5.64 


213. 


189. 


132. 


187. 


129. 


110 


127. 


0.044 


12.2 


-0.172 


6.74 


237. 


215. 


148. 


214. 


147. 


116a 


151. 


0.067 


16.8 


-0.097 


6.46 


249. 


206. 


121. 


206. 


121. 


116^3 


88.8 


0.055 


8.80 


-0.214 


8.11 


167. 


155. 


92.6 


153. 


92.0 


164 


111. 


0.038 


9.63 


-0.193 


6.88 


199. 


186. 


120. 


184. 


118. 


165C 


228.* 


0.043* 


7.70 


-0.136 


7.43 


176. 


160. 


92.7 


160. 


92.6 


211 


92.7 


0.095 


7.29 


-0.189 


7.15 


153. 


136. 


70.8 


135. 


69.0 


212 


130. 


0.062 


13.2 


-0.155 


5.64 


243. 


213. 


132. 


213. 


131. 


633A 


156.* 


0.067* 


19.9 


-0.161 


5.78 


265. 


227. 


160. 


227. 


161. 


643a 


158. 


0.086 


20.6 


-0.118 


5.78 


275. 


222. 


153. 


222. 


152. 


643/3 A 


154.* 


0.067* 


20.3 


-0.184 


7.43 


234. 


209. 


149. 


204. 


144. 


643/3B 


154.* 


0.067* 


15.1 


-0.166 


8.11 


170. 


157. 


110. 


156. 


108. 


6437 


116. 


0.054 


12.7 


-0.174 


8.53 


205. 


176. 


105. 


171. 


100. 


106 


108. 


0.075 


12.4 


-0.167 


5.64 


177. 


169. 


116. 


167. 


110. 


324 


95.2 


0.045 


8.53 


-0.201 


8.66 


160. 


147. 


89.9 


147. 


88.8 


342 


140. 


0.018 


15.1 


-0.152 


9.76 


235. 


220. 


153. 


217. 


149. 


444 


133. 


0.055 


14.0 


-0.135 


6.74 


225. 


207. 


141. 


207. 


139. 


454 


107. 


0.034 


10.6 


-0.154 


6.88 


205. 


197. 


141. 


197. 


141. 


556 


130. 


0.063 


15.4 


-0.167 


6.46 


248. 


226. 


160. 


222. 


155. 


631 


107. 


0.025 


9.63 


-0.178 


6.33 


199. 


196. 


140. 


196. 


141. 



Distances are given in kpc; velocities in km s ^. Asterisks stand for double or multiple haloes, so that their global parameters 
(with an asterisk in this Table) might not reflect a direct property of the corresponding DLO. 
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Figure 1. Upper (a): The cosine of the angle formed by the position and velocity vectors of each baryon particle within halo 
#242, versus their positions. Dots: gas particles, stars: stellar particles; open symbols: counter-rotating particles. Lower (b): 
Specific angular momentum component along j^^^ for each baryon particle of halo #242, versus their positions. Solid line: 



Figure 2. Specific angular momentum at 2 = versus the mass for the haloes (open symbols), for the inner 83 per cent of the 
gas component in the DLOs (filled symbols), and for their stellar component (starred symbols). Squares, triangles and circles 
correspond to DLOs formed in SI, Sib and Sic simulations, respectively. The solid (dotted) box encloses the region occupied 
by the spiral discs (ellipticals), as given by Fall (1983). 



Figure 3. Thick lines: the histograms for the integrated projected mass density of baryons, M=y'(Hi), for DLOs #242 and 
#643/3A (see Table Thin lines are their best fits by means of a bulge-disc decomposition. A solid-line arrow marks iJtrana 
and a dashed-line arrow marks Rmin- 



Figure 4. The bulge and disc scale lengths, iJb and R^, for double exponential fits. Squares, triangles and circles correspond 
to DLOs formed in SI, Sib and Sic simulations, respectively. Dots are from ID decompositions of surface brightness profiles 
given by CdJB96 and Courteau (1997a). The error bars in the top left corner are the typical dispersions in the i?b Eind Rd 
values for DLOs, obtained from Monte Carlo realizations of their profiles. 



Figure 5. Average values of the iJb and parameters and their dispersions, obtained from Monte Carlo realizations of 
DLOs #242 and #545 S"sht(ij) profiles (see text). Thick lines correspond to ^ = 1 (that is, Sl's'"(i?) oc S™=-==(i?)), thin 

lines to = 3 and dotted lines to = 5. 

Id Id 



Figure 6. The mass bulge-to-disc ratios B/D versus the optimal bulge shape parameter, riopt, for the DLO sample. Squares, 
triangles and circles correspond to DLOs formed in SI, Sib and Sic simulations, respectively. 
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Figure 7. Points: the rotational velocity distributions for the baryon particles of DLOs #242 of SI, #212 and #643/3A of Sib 
and #342 of Sic (see text). Solid lines: circular velocity curves (see Eqs. (^) and (^)) for the total mass distribution; dotted 
lines: the baryonic contribution Vbar('")i a^nd dashed lines: the dark halo contribution Vdm('")- 



Figure 8. The disc scale length versus Ar200 for the DLO sample. Symbols are as in Fig. ^ DLOs belonging to haloes 
that include several objects (either DLOs, spheroids or irregulars) have not been considered. 



Figure 9. The V2.2 velocities versus the disc scale lengths for the DLO sample. Squares, triangles and circles correspond to 
DLOs formed in SI, Sib and Sic simulations, respectively. Dots are data from the Courteau sample (Ha spectra, see C97b). 
Open circles are data from spirals with Hi rotation curves from CvG91 and Rhee (1996). 



Figure 10. The f-^^) versus V200 for the DLO sample. Symbols are as in Fig. 0. DLOs belonging to haloes that 

V ^cir / R2 2 

include several objects (either DLOs, spheroids or irregulars) have not been considered. The shaded area shows the 95 per cent 
confidence interval for maximum discs (see text). 



Figure 11. The rcross parameter divided by the virial radius r2oo versus the V2.2 velocities for the DLO sample. Symbols 
are as in Fig. ^. DLOs belonging to haloes that include several objects (either DLOs, spheroids or irregulars) have not been 
considered. 



Figure 12. The logarithmic slope, LS, versus the disc scale lengths (a) and the V2.2 velocities (b) for the DLO sample. 
Squares, triangles and circles correspond to DLOs formed in SI, Sib and Sic simulations, respectively. Open circles are data 
from spirals with extended Hi RCs, in the same dynamical range (see table 2 in CvG91). 
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Figure 13. The peak circular velocities yf"^^ versus the -Rpoak/Rd ratios for the DLO sample. Squares, triangles and circles 
correspond to DLOs formed in SI, Sib and Sic simulations, respectively. Dots are data from the Courteau sample of RCs in 
Ha; open circles are data from spirals with extended Hi RCs from table 4.1 in Rhee (1996). See text. 
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